Width of homoclinic zone for quadratic maps. 

V. Gelfreich f, V. Naudot * 

Mathematics Institute 
University of Warwick 
Coventry CV4 7AL, UK 

email: f V.Gelfreich@warwick.ac.uk, 
* V.Naudot@warwick.ac.uk 

May 6, 2008 



'Corresponding author 



1 



Abstract 



We study several families of planar quadratic diffeomorphisms near a 
Bogdanov-Takens bifurcation. For each family, the associate bifurcation 
diagram can be deduced from the interpolating flow. However, a zone 
of chaos confined between two lines of homoclinic bifurcation that are 
exponentially close to one-another is observed. The goal of this paper is 
to test numerically an accurate asymptotic expansion for the width of this 
chaotic zone for different families. 



Mathematics Subject Classification: 37D45, 37E30, 37G10. 



1 Introduction 

In this paper we study homoclinic bifurcations in the unfolding of a diffeomor- 
phism near a fixed point of Bogdanov-Takens type. To begin with, we consider 
a planar diffeomorphism F : ^ with the origin as a fixed point and where 

dF{0,0) ^ Id + N 

where N^O is nilpotent. The origin is said to be a fixed point of Bogdanov- 
Takens type. This latter terminology is more known for a singularity of a vector 
field X with linear part having double zero eigenvalues and a non vanishing 
nilpotent part. Since this singularity is of codimension 2,. i.e., is twice degen- 
erate, a generic unfolding will depend on two parameters say (/i, v). In the case 
of a vector field, such unfolding has been studied in [TJ [33] and for maps in 
[8l[9]. For completeness, the corresponding bifurcation diagram is revisited in 
Figure 1 on the left: a curve of homoclinic bifurcation emanates from the origin, 
below a curve of Hopf bifurcation, see [Tl] for the terminology and more details. 
For parameters located between these two curves, the corresponding dynamics 
possesses a stable limit cycle. Finally, for parameter on the ordinate {/i = 0}, a 
saddle node occurs, see also [11] for more details. 

The Bogdanov-Takens bifurcation plays an important role in dynamical sys- 
tems, for instance from the bifurcation theoretical point of view. Given any 
dynamical systems depending on a parameter, the structure of the bifurcation 
set can be often understood by the presence of several high codimension points 
which act as organising centres. Knowing the presence of (degenerate or not 
degenerate) Bogdanov Takens points initiate the searches for subordinate bifur- 
cation sets such as Hopf bifurcation sets or homoclinic bifurcation sets. In this 
paper, we consider a nondegenerate Bogdanov Takens point. 

For the map F, an unfolding theory is developed in [SKH]. It is very similar 
to the case of a flow. To be more precise, any unfolding 

of the map F (where (/i, j/) e ^^^^ p ^ p^^^^ 

can be embedded into a 

nonautonomous and periodic family of vector fields X^.u- The diffeomorphism 
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coincides with the time 1 map of that vector field, see also [33]. Using an 
averaging theorem [28j the dependence on time is removed to exponentially 
small terms. Moreover, one can show that is formally interpolated by 

an autonomous vector field see [19]. This latter can be used to study 

the bifurcations of fixed points of F^.i/. Both approaches move the difference 
between these two types of bifurcations beyond all algebraic order. 

Although all the Taylor coefficients of can be written, there is no reason 
to expect convergence of the corresponding series, since the dynamics for a 
planar diffeomorphism can be much richer than the dynamics of a planar vector 
field. In the real analytic theory, this difference is exponentially small [lOl [19] . 



Figure 1: The Bogdanov Takens bifurcation for a flow (left) and for a diffeo- 
morphism (right). 

As we said above, for diffeomorphisms, the bifurcation diagram (figure 1, on 
the right) is essentially the same. However, there is no reason to expect a sin- 
gle homoclinic curve, since a homoclinic orbit may be transverse and therefore 
persists. We observe a separatrices splitting and instead of a single homoclinic 
curve, one observes two curves J^^(/i) and i^^(/i) respectively corresponding to 
the first and the last homoclinic tangency. If a parameter (/i, ly) is (strictly) 
located in region between those two curves, then the map F^^^^ possesses trans- 
verse homoclinic trajectories. On the lower and upper boundary the homoclinic 
connexion becomes non-transverse. Understanding the width of this region is 
the main goal of this paper. 

Before going any further, we set the following preliminaries. Without lost of 
generality and up to an analytic change of coordinates one has: 



V 




xi=x + y, yi = ?/ + /p,^(x,y). 



(1) 



where 



/o.o(O) = 



dx 



(0,0) 



dy 



(0,0). 
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We shall assume that 



^(0,0)^0. (2) 
By the implicit fmiction theorem, there exists z^) such that 



dx 



(i(/i,j.),0) = 0. 



Applying a conjugacy of the form x — x + x^^^, y = y (and after removing the 
bars) amounts to writing 

fii,i^{x,y) = -boQ{fi, i^) + b2Q{fi, i^)x^ + bQi{fi, i^)y 

+ bii{^i,iy)xy + h.o.t{x,y) (3) 

h.o.t(a::, y) stands for the higher order terms in x and y. From ^ 620(0, 0)7^0. By 
a linear rescaling in the variables {x,y), we can fix 620 (m,^^) = 1- Furthermore, 
we put 611(0,0) = 7 and assume that the map 

df 

(/x,^) ^ (_/^_,(0,0),4f^(0,0)) = {-boo{fi,<^),bo,{fi,,^)) 
dy 

is a local diffeomorphism near (0, 0). From now on, we shall consider (600, 601) as 
our parameters and rename them (again) by (/x, v) i.e., write (600 1 601) = (/x, v). 

In [S] it was shown that 

^*(m) - ^-{l-2)^+0{^Jii). (4) 
In |19j the following formula is proposed: 

z/+(/i)-j.-(/i) = @^K[^Ji,^-2) + o{^?'Hog^J), (5) 

where 

K{^,, 7) = — . e-^^'/ (6) 

is referred to as the 'leading part' of the width and 0-y is an analytic invariant 
of the map i^o.o called a 'splitting constant', see pU] . 

The goal of this paper is to establish, numerically, a more accurate formula 
for the width of the homoclinic zone v'^ (/i) — v~ (/i) . The existence of asymptotic 
expansions for the width of the homoclinic zone is unknown so far. Furthermore, 
if it does exist, it is very hard to compute analytically. The difficulty here 
comes from the fact that the normal form of the map coincides with that of 
the time 1-map of a vector field. Therefore the difference between the flow 
and the map is pushed beyond any algebraic order. In the nearly integrable 
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context, a polynomial asymptotic expansion for the splitting of the separatrices 
is proposed in [29]: the author considers the perturbation of a Hamiltonian 
(elliptic) billiard. The system depends on a perturbation parameter e > 0, a 
hyperbolicity parameter h > and admits four separatrices, which break up 
when h > 0. In this special case, the author proposed an asymptotic expansion 
for the area of the main lobes of the resulting turnstile that takes the form of a 
power series (with even terms) in e. See [211 [HI [M] for more references on the 
computation of separatrices splitting. 

In this paper our approach is somehow experimental. We study examples 
and present strong numerical evidence for the following expansion of the width 
of the homoclinic zone 

^+(M)-^-(/i)xif(^,7-2)^ ^ Sfe,,//4logV (7) 

fc>Oo<j<[|] 

k 

where [— ] stands for the integer part of fc/2 and K is given by (O. 
Remarks: 

(a) Observe that ([7]) is a double series with logarithmic terms and numeri- 
cally, for such an expansion, we do not know any efficient techniques to compute 
the corresponding coefficients with a large precision. However, our numerical ex- 
periments showed that log{i>'^ (fj.) — v~ (fj,)) has a simpler asymptotics expansion, 
than J/+ (/i) — (/i) itself. More precisely we have 

log f iy+ (/i) - ly- (/i) j X log f K{fi, 7 " 2) j + ^ rukfi^ + log /i ^ nfe^^ . (8) 

One easily checks that formula ([7]) follows from ([8]) and that the Cfej 's depend 
on the TOfc's and the rtfe's. Note that the asymptotic series ([8]) does not involve a 
double summation and therefore the corresponding coefficients can be computed 
with a much higher precision. 

(b) Logarithmic terms may vanish, this occurs for instance in the case of the 
Henon map, see next section for more details. 

(c) From the numerical data, we are able to guess a simple analytic expression 
for the first logarithmic term in ([8|). More precisely we have 
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which is valid for all families studied in this article. The paper is organised as 
follows. We shall consider three different families that satisfy, (up to appropri- 
ate smooth changes of coordinates) the setting above with different non linear 
terms. As a result of our experiments, for each family we shall state the asymp- 
totics for the width of the homoclinic zone, confirming formula Looking 
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for the width of the homochnic zone amounts to fixing one parameter, say ^ 
in the unfolding ([3]), and find the values of the second parameter, say for 
which the system admits a first and a last homoclinic tangency. We say that 
is the main 'parameter' and v is the 'slave' parameter. In section 3 we briefly 
present the strategy to follow. The rest of the section is devoted to the com- 
putation of the invariant (stable and unstable) manifolds at the saddle point. 
The splitting function which is a key ingredient of the techniques is presented. 
Indeed, primary homoclinic orbits are in one to one correspondence with zeroes 
of the splitting function. Therefore, the first and the last homoclinic tangencies 
will correspond to double zeroes of the splitting functions. Moreover, the split- 
ting function is periodic, with exponentially decreasing harmonics and is well 
approximated by the splitting determinant. With a good precision, computing 
the width of the zone amounts to the computation of the first two harmonics of 
the splitting function and their dependence with respect to the slave parameter 
(the main parameter being fixed). For each family, we compute the width of 
the homoclinic zone for several hundreds values of the main parameter and 
collect the results in a set of renormalised data. In the next step, the coefficients 
in ([8]) (considered as an ansatz) are extracted by interpolation techniques. The 
remaining part of the paper is devoted to the verification of the validity of our 
results. More precisely, we test the ansatz ^ and we find how precise our 
data for the width of the homoclinic zone should be in order to produce reliable 
results for the coefficients of the asymptotic expansion. Finally, the constant 
coefficient of the expansion should coincide with the splitting constant [T9j : fol- 
lowing the procedure developed in [20], we compare these constants with the 
constant coefficients of the expansions. 

2 Main results 

Before presenting our main results, we first introduce the following notions. 



2.1 Asymptotic sequences and expansions 



Let Eq > be given and let 

S = {fo: fl, ■ ■ ■ , fni ■ ■ ■} 

where /o = 1 and for each integer i > 0, fi : (0, eo) ^ M is a smooth positive 
function such that 

or in other words /i+i(x) = o{Ji{x)). Such a family S is called an asymptotic 
sequence. In this paper we shall consider the following asymptotic sequences 

V = {l,a:,a;2...,x",...} (9) 
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that is fi{x) — and the Dulac asymptotic sequence [57]: 

P = {1, X, logo;, a;^ a;^ a;"* logo;, a;^ . . . , x^" log(a;), a;^", x^^+i, . . .} 
that is for all integer n > 

/3„(X) = X^", hn+l{x) = X^^+\ hn+2{x) = X^^+^ \og{x) . (10) 

Let (j) : (0,eo) — > M be a smooth fmiction. We say that 

m X ^a„/„(x) (11) 

is an asymptotic expansion of (/> at (where the {/n}neN is an asymptotic 
sequence and all a„'s are real) if for all integer n, 

n 

cj,{x) - 0{">(x) = 0(/„+i(a;)), where ^^(a;) = ^a,/,(x). 

When looking at expansions of the form (llip no convergence is implied and 
often the a^'s are Gevrey-1, i.e., 

3 M > 0, r > 0, such that Vfc > 0, \ak\ < Mk\/r^. (12) 

2.2 Quadratic family 

Our first example is the Quadratic map 

{x,y) ^ {x + y,y + x'^ - n + -yxy + vy) 

Observe that Q mimics the unfolding ^ i.e., that takes the form of F^^i, and 
ignores the higher order terms. We normalise the width of the homoclinic zone 
associated to the Quadratic family by defining 

Q (,,\ - '^'^(m) - '^"(m) 

Within precision of our computations we observe 

logS.ifi) X ^Affe(7)//^ + logA*^iVfc(7)/i'/', (14) 

fc>0 k>l 

where Mk{'y) and Nk{j) are real coefficients which depend on the parameter 7. 
Comparing with ([5]), we see that 

exp(Mo(7)) = &j 
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is the splitting constant associated with Qo,o,7- Moreover, as we announced in 
the previous section, we have 



/ 6(7-2 )-^ 



For each value of 7, the M^'s and A^^'s can be computed with a very high 
precision, see Table 1 for illustration. Formula p4|) is verified for the 76 first 
coeficients: Mk, k — 0, . . . , 50 and Ni, £ = 1, . . . , 25. Although the precision 
decreases almost linearly as k and £ increase, the 76 first coefficients can be 
computed with 60 correct digits. To compute these first coeficients, we need to 
compute the width of the homoclinic zone with at least 200 correct digits, see 
section 4.2 for more details. 

Even if we can propose an analytic expression for iVi(7), we have not been 
able to guess analytic expressions for the other coefficients Nk and Mk- 

2.3 Bogdanov family 

Our second example is the Bogdanov map [3l [H [6] . 

{x, y) I— > (x + y + + ^xy + ax + by ,y + + jxy + ax + by). 

The Bogdanov map, see for example [51|7j, is the Euler map of a two-dimensional 
system of ordinary differential equations. In [3] Arrowsmith studied the bifurca- 
tions and basins of attraction and showed the existence of mode locking, Arnold 
tongues, and chaos, see also [1] for more details. 

For this map the saddle point is located at the origin. This map can be 
transformed to the form ([3]). Indeed, let 

u — X — a/2, V — y + {x — a/2)^ + 7(2: — a/2)y + a{x — a/2) -f by. 

We retrieve the map (2) and higher order terms ^ by putting 

u^a + b-{^ + 2)^, 7 = 7 -f 2, ^ = aV4, 



and 



ft.,y = {x + y)^ ~ ^i + jy'^. 



The parameter a is chosen to be the main parameter and b the slave parameter. 
From ([4]), the Bogdanov map admits a homoclinic zone near the line 

6±(a) = ^aj + 0{a^/^). 

The normalised width takes the form 

b+{a)-b-{a) 



S^{a) 



^(aV4,7) ■ 
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coef. 


scale 


value 


coef. 


scale 


value 




1 


61.26721889 




1 


-13.35083105 


Ai 




-29.82701974 


Ml 


1 //I 


-35.34533603 


Bi 


a log a 


-6.612244898 


Ni 


^1/2 log^ 


-9.183673469 


A2 


a 


5.824479250 


M2 




-25.71572403 


A3 




17.41183781 


M3 




60.69366755 


B2 


log a 


5.649967276 


N2 


/ilog/i 


-41.92449575 


Ai 


a' 


-0.2874798361 


M4 


M 


-215.4221683 


A^ 




-22.04012159 


Ms 


r / .1 


-45.92851439 


B3 


a"" log a 


-6.966574583 




^3/2 log^ 


-242.5333437 


Ae 


of- 


-6.250578833 


Me 




-960.8699623 


A7 




39.27382902 


My 


TTT 


755.3601690 


B4 


log a 


10.92891913 


iV4 




-1587.303140 


As 




19.31687979 


AIs 




-3308.441120 


A9 




-82.17477248 


Mg 




1090.837521 


Br, 


log a 


-20.01663759 


N5 




-11017.80445 


Alo 




-50.35178499 


Mio 




-134120.3771 


An 




186.9039750 


Mil 




22519.75418 


Be 


a'' log a 


40.63376347 


Ne 


log A* 


-79363.78673 


Aw 




128.7996196 


Mio 




904656.6104 


All 




-444.7385574 


Mil 


^13/4 


87833.05069 



Table 1: The 20 first coefficients of the asymptotic expansion for the Bogdanov 
map (left, 7 = 8) and the Quadratic map (right, 7 = —3). All the given digits 
are correct. 

Similarly to the Quadratic family, our experiments showed that log satisfies 
the following asymptotics: 

log^^(a) X ^A-(7)a'=/'+loga^Sfc(7)a^ (15) 

k>0 k>l 

where ^^(7) and Bk{'y) are real coefficients which depend on the parameter 7. 
Comparing with (O, we see that 

exp(Ao(7)) = 0^ 

is the splitting constant associated with Bo.o,7- Moreover, we observe numeri- 
cally that Bi(7) = -(67/7)^ 

In Table 1, we provide typical results for our computation for the Quadratic 
and Bogdanov maps. Although the first 20 coefficients do not show a tendency 
to grow rapidly, we conjecture the series and (fT5)) diverge and belong to 
the Gevrey-1 class (IT^ . compare with pT) . 



9 



2.4 Henon map 

The last example to be considered in this paper is the Henon map [52] defined 

by 

H = H-j; : {u,v)^ {ui,vi) 

where 

Ml = V, vi = av^ — hu + 1. 

See |23) for recent results concerning this family. The Henon map has a fixed 
point of Bogdanov Takens type at a = 6 = 1. We chose a as the main parameter 
and b as the slave parameter. We note that the Henon map is conjugate to the 
Bogdanov family in the special case of 7 = 0. The conjugacy is given by the 
following change of coordinates and parameters 

u = x, V = X + y + + ax + hy, 6 = 6 + 1, a = (1 + b/2f - 

We also observe that the Henon map can be transformed to the form ^ with 
the non linear term of the form ([3]) by putting 

I, 6+1, I, 6+1, 1 

u = -(a; + ^5—), v^-[x + — — ) + -y. 
a I a 2 a 

In the new system of coordinates, the Henon map takes the form (|3]) with 

^ (x + yf - + vy, 

where 

The Henon map admits a homoclinic zone near the line 

b^{h) = 1, 5 > 1. 

In the case of the Henon map we define the normalized width of the zone by 

~ b+{h)-b-{a) 
= K{l-~a,0) ■ 

Our numerical experiments show that S has the following asymptotic expansion: 

S{a) = ^ife(l-a)'=/4. (16) 

fc>0 

Unlike the case of the Bogdanov map with 77^ 2 (i.e., 77^ 0), the asymptotic 
expansion does not contain logarithmic terms. We expect this property to be 
closely related to the fact that the Henon map contains a one parametric sub- 
family of area preserving maps. In general, even when 7 = 2, there is no reason 
to expect the logaritmic terms to vanish for a map F^,^. 
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coef. 


scale 


Hcnon map 


Ao 


1 


2.4744255935532510538408* 10 


Ai 


|a-l|i/4 


-2.878113364919828141704* 10^ 


A2 


|a-i|i/2 


1.8211174314566012763528* 10® 


A3 


|a-lp/4 


-412552.07921345800366019 




|a-l| 


-309961.28583121907079391 


^5 


|a- i|5/4 


257055.93487794037812901 


As 


|a-lp/2 


-56830.201956139947433580 


A7 




- 12386.990577003086404843 


^8 


|a-ip 


-11792.964908478734939516 


^9 


|5-l|9/4 


18742.189161591275288347 


Aw 


|a-l|5/2 


-4774.6727458595190485600 


All 


|a-l|"/4 


-2822.9663193640187675835 


A12 


|a-i|3 


3276.6438736125169964394 


Aw 


|5-l|"/4 


-1910.5466958542171966392 


Au 


ia-ir/2 


7704.6605615546853854041 


Aw 


|a-i|i5/4 


-7827.0351891507566506398 


Aw 


|a-i|4 


13919.102717097324631620 


An 


|a-i|i^/4 


-11932.139780641352182621 


Aw 


|a-l|9/2 


22120.721696311178434645 



Table 2: The 19 first coefficients in (fT5|) . All the given digits are correct, 
also conjecture that the series ([TSl) belongs to the Gevrey-1 class. 
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3 Computing the width of the homocUnic zone 



In this section, our approach concerns the Quadratic family Qp,,y.-y. The other 
famihes (Bogdanov and Henon) are treated in a similar way. From now on, we 
do not mention the (/x, 7) dependences when it is not necessary, but we may 
emphasise that dependence when it is needed. 

3.1 Strategy 

i) We assume an ansatz and in particular the one given in formula ([8]); 

ii) Compute n (several hundreds) values of the width for values of /i^/** € 
[c, d], where < c < d are close to (typically c w 5/1000, d w 1/100). It 
is convenient to work with the so called 'normalised width of homoclinic 
zone' defined by 

X(A.,7-2) 

where K is defined by The result is collected in a set of data of the 
form 

H ^ {{fif ^og{S^{^i^)), c<fi^<d, i = l,...,n}. (17) 

iii) Take ^ € N such that 3^/2 + 1 < n and i >> 1 even. Then we compute 
the coefficients Mk, k = 0, . . . ,£ and Nk, k = 1, . - . , £/2, of the truncated 
expansion 

G{3^/2>(A.) = ^M,(7)//^ + log/if: A^.(7)//^ 

fe=0 k=l 

to interpolate the set H, i.e., for all integer i = 1, . . . , 3£/2 + 1, we have 

k=0 fe=l 

See subsection 3.10 for more details. 



Remarks: 

• For the Bogdanov family, the set of data for the normalised width is de- 
noted by 

ii = {{aj ,\og{Sj{ai))), c < Qi < d, i = l,...,n}, (18) 

where 



~ b+{a,)-b (g,) 

bj[ai) = — „ _ — and < c < rf. 
i^(a|/4,7) 
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• For the Henon family, the set of data for the normahsed width is denoted 

by 

Z = {{\\-~a^\,S(a,)), c<\^-~a^\<d, ^ = (19) 

where 

^ ^» = i^ii -An and < c < d 
3.2 Invariant manifolds 

We now compute the stable and unstable manifold at the saddle point. In what 
follows, our description concerns the Quadratic map Q but similar computations 
are done for the Bogdanov map and the Henon map. 
From (|13p . the map Q has two fixed points 

= (VTi^O), and = (-^71,0). 

is a focus and is a saddle and will be the point of interest. The eigenvalues 
of (iQ(S^) are given by 

Ai = ^(^2 + i. + 7Vm- ^(7VM + i/)2 + 8VM^, 

A2 - ^(^2 + i. + 7Vm+.^(7Vm + //)2 + 8V7Z^. 

For /X > sufficiently small it is clear that Ai < 1 < A2. At the saddle S^, the 
Taylor expansion of the local stable manifold W(^^ and that of the local unstable 
manifold Wj"^ are computed as follows. Denote by 

y 00 00 

: (M, 0) ^ (K^ S^), z ^ $,(z) = f + ^ ^kz\Y. ^^^^ 

^ fc=i fc=i 

r 00 C30 ■ 

^ k=\ k=\ 

the parameterisations which respectively satisfy 

$,(Aiz) = Qo$,(z) and $„(A2z) = Q o $„(z) (20) 

for all z near 0. Substituting the series into (PO)) and collecting terms of the 
same order in z we get 



Ej=o '/'i'/^'fc-i + 7 Ej=o '/'i'^fe-j + ^^k = A^-0fc 

Pk + fk = A^pfc, fc>l 
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Since A2 > 1 and Q is entire, from (^0)) we easily deduce that the radius of 
convergence of the series defined in ((22|) is infinite. Denote by g the radius of 
convergence of the series defined in (pij) . We fix N^^^ e N. Since we are after 
a single branch of the stable manifold we write 

where 

^ fe = l 

and where Q < 5s < Q- We proceed in the same way for the local unstable 
manifold, i.e.. 



*-..iv„ax(^) = Vm+ E /fe^"' E Pfe^' (23) 



where 



k=l k=l 

and where < 5^ << I. The local invariant manifolds are computed with the 
following precision: 

\\^.,N,..Jz) - ^siz)\\ = 0(2^—), ||$u,JV„ax(^) - '^u{z)\\ = ). 

In particular we have 

||<i>„(A2z) - Q o <i>„,jv.„ax(^)ll = ). (24) 

Since we need to study the map when homoclinic orbits are present, we need 
a good estimate of the global unstable manifold. Recall that $„ is entire and 
therefore both components defined in ([23| converge for all z as N^^x 00. 
However, for large z, the computation of the unstable manifold requires too 
many coefficients and therefore p3| is not very convenient. We then proceed as 
follows. Let Pq — $«(zo) & and choose mp such that 

Then, for any fixed mo, we have 

Po = lim Q™°o$„,jv„,ax(^i) 

and if zi << 1 the convergence is fast. Therefore, by putting 

- ^^„ax,m = {Q" ° *u,A^»ax(A^"2), < z < zo}, m > mo 
we get an accurate estimation of the global unstable manifold. 
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3.3 Jacobian and Wronskian functions 

Before introducing the sphtting function which wiU play a key role in the paper, 
we need to introduce two additional functions. We first define 

J : P ^ C, z detdQ($s(z)) 

as the Jacobian of the map Q along the stable manifold 

A straightforward computation gives 

J{z) = 1 + + (7 - 2)$,,,(z) - 7$.,y(^)- (25) 
In terms of series, from we get 



J(z) = ^ Jfcz*=, where (26) 

fc=0 

Jo = 1 + + (7 - 2)^, and Vfc > 0, Jfc = (7 - 2)0fe - -fi^k- 
The Wronskian function (along the local stable manifold) 

n : I? ^ R, z n{z) 

satisfies 

n{Xiz) = J{z)n{z). (27) 
We put ilo = 1 and look for a solution of (I27p of the form 

17(^) = 2l°g-^o/logAi A ^g^^^fcY ^28) 
^ k=l ^ 

With (HZl), (Uni), and it follows that 



1 ( ""^ \ 



Both series (|26p and (pS)) are convergent. The functions J and f2 will be ap- 
proximated by 

Jjv„,^Jz) = 2 Jfez*^ and f^w„„x(^) = z'°S'^°/i°g^M 1 + ^ ilkz'' 

k=0 ^ fc=l 

respectively. In this way, we have 

\^N^.A^iz) ~ JN^Jz)nN^^^{z)\ = ). (29) 
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3.4 Splitting function and flow box theorem 



In this section, we introduce the key part of our techniques. Recall that in 
our investigation for the width of the homoclinic zone, we fix the value of the 
main parameter and look for values i/+ and v~ of the slave parameter that 
correspond, respectively, to the first and the last homoclinic tangency. In order 
to find a homoclinic point we need to adjust the slave parameter in such a 
way that two curves on the plane have an intersection. Finding a homoclinic 
tangency requires additional adjustments to make this intersection degenerate. 
This problem is much easier in the discrete flow box coordinates, in which the 
stable curve coincides with the horizontal axis and the unstable one is a graph 
of a periodic function. A further simplification will be achieved by observing 
that this periodic function is very close to a trigonometric polynomial of the 
first order. The splitting function ~ we shall introduce now is such that 
the first and the last tangency correspond to double zeroes of 0^^;^+ and 6^, ir- 
respectively. Our investigation amounts then to finding values z/"*" and v~ such 
that 0^,;y+ and 0^,,^- possess double zeroes. 

In this section, we present the splitting function 0^,,^ for the Quadratic map, 
in the case of the Bogdanov map, the splitting function is denoted by 6a. 6. In 
what follows, we assume that the parameter (^, v) is such that the map Q 
possesses a homoclinic orbit, i.e., the unstable manifold intersects the local 
stable manifold at a point = (70 = 'I's(-Zs)- Then we fix a neighbourhood 

U of the point go- We parametrise W(^^ near go by 



where Iq = (—1,1) and near go by 



r„ : lo <I>„(2:„ • A*). 



Now we state the following (flow box) lemma [13] . 



Lemma 1 There exists i?o > and an analytic diffeomorphism 



* : {-Eo,Eo) X lo 
{E,t) 



2 



I — > 



such that the following hold 



i) 'i>{E,t + l) = Qo^{E,t), 



n) ^-(0, 0) = go, *(0, t) G W{^^ for teh, 
iii) the Jacobian matrix 




(30) 



is such that the second column of (i^(0,<) is Ts = dTs{t)/dt, 
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iv) the map fl{E,t) = det d'^{E,t) satisfies 0(0, = fl{zs ■ X\); 

The splitting function, denoted by Ofj,^^{t), is the first component of 

^-^oTuit)-^-^oT,{t). 
Applying Taylor theorem at the stable manifold, we get 

+ o(^\\T^{t)-TS)r)- (31) 

The following properties hold: 

[-] Let < S < TT. The map Qn,v has an analytic continuation onto the 
rectangle: 

B = {i e C I t = t' + it", i' elo, \t"\ < g}, \g\ < (tt - 5)/| log Ai|. (32) 
The function G^,,^ is periodic so we can expand it into Fourier series: 

oc 
j=-oo 

As usual, the Fourier coefficients are defined by an integral: 

Pk{n,u)= [ @^At)e-'^'^'^*dt, for each fceN. 
Jo 

Let < e < (tt - log(Ai)|. Since the integral of 6^,i.(i)e~^*'"'* over the 
boundary of the rectangle {{£ + Of') | < t' < 1, < t" < £>} vanishes, we 
conclude 

e^,,(t)e-2''='^*di = e-^'^^e e^_^(i + iQ)e-^^^'''dt. (33) 

JO ^0 



Consequently 



|Pfc(M,J^)| <sup|e^,,(t + ie)|-e-2|'=l-^ (34) 



i.e., the harmonics of decrease exponentially. The function 0^^^ can be 
well approximated by the sum of zero and first order harmonics: 

e^,,(t) = P_i(M,!^)e-2'"* + Po(M,i^)+Pi(M,i^)e''"* + 02(0 (35) 

or equivalently, Qix,v is well approximated by a trigonometric polynomial func- 
tion 

®nAt) = Po(m, + 2|P_i(m, v)\ cos(27rt + arg(P_i(/i, v))) + 02{t) (36) 
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where 



SUPIO2IW = o(sup|e^,,(t)n. 

iGlo telo 



(37) 



[-] Since d4'-i(^'(0, t)) = {d'^{0,t))-\ we have 



Furthermore, 

*-'(r"(t)) = (^ii;„(i),T«(i)), *-i(r(0) = (^i?.(i),T.(t)) = (o,<), 

with dSOj) and jST]) it follows that 

Q^.A^) = E^{t)~E,{t) = _l_det(^ |r,(t) , r„(t)-r,(t)) 

+ O(||r„(t)-r,(0ir). (38) 

Thus, we obtain a formula suitable for computation of the splitting function in 
terms of the parametrization of the stable and unstable manifold: 

e;.,.(t) = e^,.(t) + /i^,.(i) (39) 

where 

e.,.(0 = fX^'^Kl"^^^'^ r.(t)-r,(.)) (40) 

is the splitting determinant and 

IV,(t)| = o(sup|e^.,(t)n. (41) 

telo 

Note that even if the invariant manifolds and the Wronskian are computed with 
a very high precision, the function 6,j,iy(t) is only evaluated with a relative error 
of order (^(suptgi,, |e^,^|). 



3.5 Approaching a primary homoclinic orbit 

In order to compute the width of the homoclinic zone, we first find a value v = v 
where the map possesses a primary homoclinic orbit. Near v = Lemma 1 
will then be applied and the splitting determinant 6^^,,^ will be computed. We 
proceed as follows: we fix < Zs < and a section E transverse to the local 
stable manifold at -p^ = (i>s(zs). We parametrise E as follows 

S = {pi. + (0, y), -yo<y < yo} 
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where < yo << 1. For each value of the main parameter, we consider the slave 
parameter being close to I'o — (6(7 — 2)/7)y^and compute a point G W^HT, 
which is the 'first intersection' of with the section. In order to increase the 
speed of computations we use Newton's method to solve the equation r„(i) £ S. 
After that we adjust v in such a way that q^, — p^. We do not know an easy 
way to evaluate the derivative of with respect to therefore we cannot 
apply Newton's method. However, we replace the derivative by a finite difference 
approximation and use the so called 'secant' method. In other words we consider 
the limit of the following sequence: 

Vn+l — t^n H 

where qv ^ Vv + (0, Dv) and where < (5 << 1. Denote by 

V ~ lim Vn- 

n — >co 

Since pp — qp, the point (/i, v) belongs to the homoclinic zone. 

Our next step is with the computation of the width v^{ij) — (a*) for the 
given value of /x. The zeroes (and double zeroes) of 6^^^ are in one to one corre- 
spondence with primary homoclinic orbits (and homoclinic tangencies) for the 
corresponding map, see [15|, I19j for more details. We then replace the problem 
of finding homoclinic points and homoclinic tangencies by finding double zeroes 
of the splitting function 6^,1/. 

3.6 First and last tangency 

The most natural way to compute the width of homoclinic zone is to estimate 
both ~ '^^(/^) and = v^{^). Write 

Qf^At) = Po{p,'y) + et.At)- (42) 

At the first tangency, (i^ — v^) the graph of the splitting function is located 
below the t axis and Q^^v- admits a double zero. Therefore there exists G Iq 
such that 

0,.,.-(i") = supe,,,,-(t) = = Po(Ai,t^")+supe^.,-(t). (43) 

tGio teio 

At the last tangency, {v — v^) the graph of the splitting function is located 
above the t axis and 0p,jy+ admits a double zero. Therefore there exists G Iq 
such that 

e^,.+ (i+) = inf e^,,+ (t) = = Po(a*, y+) + inf e^,.+ (i)- (44) 

If we neglect O2 in ([55]) , ([^5)1 and (HH) are equivalent to 
Po(m,'^+) - 2\P^^(pL,y+)\ = 0, 

(45) 

Po(M,i^") + 2\P^i{fi,u-)\ = 0. 
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In this way the problem of finding the first and the last tangencies, is replaced 
by scalar equations in one variable each. Therefore, instead of looking for in- 
tersections between Wj"^ and and their tangencies, we save a lot of time by 
simply solving a scalar equation. Observe that for v near for all t g Iq we 
have 



Po(a*, = 0(|Pi(m, '^)), sup 16^^,(01 - 0(|Pi(m, i^))- 



tela 



(46) 



From (|55|) we need only 4 points per-period to evaluate Pq and P±i. Concretely 
we write 

Po(/i,:/)«Ro(A*,J^) = i(e^,.(0) + 9^^,(1/2)) 
F.l{^l,ly)^■R-l{^l,l^) = i(e^,,(0)- 9^^,(1/2) 

+ *(e^^.(i/4)-e^,,(-i/4))) 
Pi(/i,i/)«Ri(/i,j.) = i(e^,,(o)-e^,,(i/2) 

- *(e^,.(i/4)-e^,,(-i/4))). 

From (j35[) and (|39p . the approximation here means 



(47) 



max{|Ro(Ai,i^) -Po(A*,i^)UR±i(Ai,i^) -P±i(/i,'^)|} =0(sup|e,,.,(t)n. (48) 

teio 



Moreover, with we have 

|R±i(M,i')-P±i(M,'^)l = 0(|R±i(M,^^)|'). 

We then solve 

Ro(Ai,i>+) - 2|R_i(^,z>+)| = 0, 
Ro(A*,i>-) + 2|R_i(/x,:>-)| = 0. 
From (gSl), mi), (gHl) and ([501), we have 

e,,,,- (t-) - o(r2_i(m, i>-)), e^,,+(t+) - o(r2_i(a., ;>+)). 

By the Mean Value Theorem, we have 



O 



(49) 



(50) 



(51) 



(52) 



This approach gives a good estimation of the locus of the homoclinic zone 
and therefore of the corresponding width, but requires the computation of both 
v'^ and v~ with a very high precision. To be more precise, assume we want 
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to compute the width of the homocHnic zone for a given value of the main 
parameter with N correct digits, while the width of the zone (roughly estimated 
with formula (6)) satisfies 

10^^+1 < i/+ - i.- < 10^% (53) 

where Nz >> 1. Thus we need to compute both and i/^ with N^ + N correct 
digits. We observe (numerically) that 

also compare with ([55]) below. Therefore with ([5^ and (|54p . i>+ — gives 
an estimation of the width with a relative error of the same order. In partic- 
ular, this means that we cannot choose N bigger than N^. With this method, 
thanks to ([49]), the estimations of Po(m, and of P_i(^, i^) are obtained with 
a relative error of the same order as |P_i(^, This requires the computation 
of the splitting determinant with the same relative precision. When the main 
parameter tends to 0, since the eigenvalues Ai and A2 tend to 1, the number of 
iterations (i.e., mo) and the number of terms in (j23p . (i.e., iVmax) required to 
compute the unstable manifold need to be chosen bigger and bigger. Moreover, 
in order to guarantee (|^ . we need to have Po(m7 = C(P-i(a*) '^)): i-e-, (^5]) . 
which requires that the local stable and the unstable manifold are close to one 
another and more precisely 

||r„(i)-rn(0)ll -o(i^(M,7-2)). (55) 

As a conculsion, when the main parameter tends to 0, this approach becomes 
more and more delicate. 

In what follows, we propose another approach which does not require the 
computation of Po(/i, 1^), still requires a first value of v — P such that (|46p and 
gives an estimation of the width with the same precision. 

3.7 'Real' approach 

From (gSl) and (gi]) we have 

Po(M,i^+)-Po(A*,i^") = -inf e,,.,+ (t)+supe,,,,-(i). (56) 

*eio teio 

Furthermore, from the Mean Value Theorem, there exists i>~ < 1/2 < such 
that 

Po(m,^^+)-Po(a^,^^-) = ^l=..-{>^+-'^~)- (57) 



Thus we get 



suPtgio ®tJ...-y- (t) - inftgio 6^^^+ (t) 



(58) 
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We observe (numerically) that 9^,^ does not change much with respect to v. 
More precisely for all v~ < < i'^ , i^^ < 1^4 < and for all t G Iq , 

|e^,.,(t)-e^,,3(t]| ^ o{\p_^\{^l,D)). (59) 



Thus, with 136]) and (gB]) we have 

supe^,,-(t)- inf e^,,+ (t) = 4\p^^\{^l,9) + oiFl,{^l,D)). (60) 

teio ■ *eio 
Furthermore with (1421) we have 



.=uAt) = ^l.=..W-^U=..W (61) 

(t) + 0(P_i(M,i^)). 



au 



With dSHl) and ([Ml) we have 



^|.=.,(i) = ^|.=..(t) +0(e^,.(t)). (62) 



We then write 



39) 

0^,.At)-Q^..s{t) = ^T^(i)|.=.3 •(^4-^^3)+0((^^4-i^3)'), (63) 

oy 



and therefore 



//4 — 1^3 dv 



{t%=u.+0{{v+ -V-)). (64) 



We observe (numerically) that the left hand side of ((64)) stays away from as 
the main parameter tends to 0, more precisely there exists > such that for 
all /i > 0, 1^3, near V and for all t £ lo, 



V4 - V3 



I > ^^0. (65) 



With ([46]) . (j6ip and ((64)) and by choosing 1/3 and sufficiently close to one 
another, we have 

OJ/ 1^4 — V3 

Therefore, with ([58]), ([601), dM]) and ([661) we can write 

V+-V- = Z{fi)+OiZ^{fi)) (67) 
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where 

Z[,) . 4|P_,|(M,P)(:.3-^4)^ 

where io S lo- 

Thanks to (|47l) we obtain the foUowing estimation for the width of the homo- 
cUnic zone 

= 4|R-i(,.,P)|(.3-^4) _ 

With (gg), (EZl) and dnSD it follows that 

{v+ ~ I.-) - Zriti) = (70) 

This 'real' approach gives a good estimation of the width of the homoclinic 
zone with the same precision as before in (|52p . Moreover, it requires only the 
computation of P_i and that of Qi_i,^{t) for two different values of v. 
However, we still need to find a value of v — v such that (l46ll holds. 



In what follows we present another way to compute the width: in the new 
approach, does not need to return near Tg as close as in (j55p . In this way, 
we will be able to compute the splitting determinant with less precision. This 
alternative approach consists of looking at the splitting function for complex 
value of t. 



3.8 'Complex' approach 

Now we present another way to compute the first harmonic, with less precision 
than in the 'real' approach case, but with less effort. Recall that formulae (5) 
and (6) already give the following estimate 

v+ -y- =0{K{^i,^-2)). (71) 

Moreover, with ([55)1 and (|7T|) gives us a rough estimate of |P±i|, i.e., we 

have \P\±^{^l,v)=0{K{^l,^-2)). 

Take Q «6 < g and Aq = K{^i,^ ~ 2)e^'^^ such that K{fi,"f - 2) << Aq. 
Assume that we have found a value oi v ^ such that 

i^(Ai,7-2) « supe^,,„(i) < Aq. (72) 

teio 

Observe that looking for such a value of = i^o requires less effort than searching 
for V where supjgj^ O^^p = 0{K{fi,^ — 2)). In particular, we only need to 
compute the splitting function with a relative error of order sup^gj^ 0p^^„(t). 
With (|35p . there exists s € Iq such that 

sup e^.,„(i) ^ Po(m, i'o) + P-i(a^, i^o)e-^'^' + Pi(m, ^^0)6^"^ + 02{s). 

teio 
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Since |P±i|(/i, t'o) = C'(i^(^, 7 — 2)), we then conclude that Po(^, t'o) = 0{Ao). 
Since Aq >> K{^, 7—2), we are not able to compute precisely the first harmonic 
P_i(/z, I'o), with the real approach. However, instead of considering t e Iq as 
real, we now consider t in the complex interval [5i,6i + 1]. Recall that the 
Fourier coefficients of Q^^^o (t) are 

Jo Jo 
Since 6^,i/o is periodic and analytic in B defined in (31), we have 

P_i(M,i.o) = / e2"*e^,,„(i)dt. (73) 

JiS 

With (123 we have 

e'"*e^,,„(t) = P_i(m, lyo) + e'"*Po(Ai, ^^o) + e'^'P^ifi, v,,) + e''^''02{t), (74) 
where 

0^{t) = 0(sup|e2,^(z<5 + <')!)• (75) 

t'eio 



With jMl) we have 

P±i(m,^^o) = 0(e-2"^). 
Therefore, since Po(m, t'o) = ©(Aq), with ((35| and ((75l) . we have 

sup |e^,,J(z5 + t') = 0(e2-^-2-f) = 0(|P_i|e2-^), (76) 

t'eio 

and further we have 

Pi(M,t^o)e^'"* = 0(e-2-(e+2A-))^ (77) 
|e2-'*O2(0l = 0(|e+2'^**|P_ip(^,i/o)e'*"*)|) =0(e-2-(2e-5)). 

We distinguish two cases 

Case 1: 5 > —. In this case, 2g — 5 < g + 2S and from ([77)1 we have 
3 

|Pi(M,i.o)e^^"*| « |e2-**02(t)|. 

Using ((39)) we write 

e'""e^,,„(t) =e2-*e,,,,„(t)+e2-*/i^,,„(t) ==^(t)+£i(t) (78) 

where A{t) = P_i(m, i^o) + e2"*Po(Ai, 1^0) and with 

Slit) = ©(e-^'^^^e-i))^ 
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Observe that for i G Z(5 + 1] 
In this case 



45+1 1 

A{t)dt = - ( A{i6) + AiiS + 1/2) ) . (79) 



But with, ([78|) we have 



1(5+1 pi5+l 

e^"'e^,.oit)dt^ / Ait)dt + Oie-^''^^e-S)y (80) 

f iS J iS 

Finally from ([73]), (Ull), (Ull) and jHO]) we get 

P_i(^,i/o) = C_i(^, z^o) + 
where C_i(Ai,i^o) = ^e-^-* (^e^,,„(i<5) - e^,,„(z<5 + 1/2) 

= 0(e-2-(2e-*)). (81) 

Case 2: i5 < -f- In this case, from ([77)) we have 

therefore we cannot neglect the term Pie"''^** from the integration in ([75)1 . Thus 
we write 



e'^"*e^..oW = P-i(M,'^o) + e2^"Po(//,i^o) + Pi(A*,i^o)e""*+f(i) 
= Ait) + Sit) 

where with ([77)1 

i(i) = P_i(/., J^o) + e'"*Po(/i, i^o) + Pie'""'', £{t) = 0(e-2-(2.-5)). 
In this case 

iS+l ^ / 

Ait)dt = - \Aii5) + Aii5 + 1/2) + AiiS + 1/4) + Aii5 + 3/4) 
and we get 

P_i(/i,t'o) = C_i(^, I'o) + ?'2, where 

c_i(/i,i.o) = e-'"*i(e^,.o(o)-e^,.„(i/2) 
- z(e^,.„(i/4)-e^^.„(-i/4))) 

|r2| = 0(e-2-(2e-5)) ^O(c^j(^,^.o)e2^^). (82) 
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When 5 > g/3, the estimation given in ([5T|) requires the computation of the 
sphtting determinant Q^^i^git) at two different values of t only. However, when 
S < e/3, requires four different values of t. The computation in the first 
case is faster, but since 6 is bigger, we loose some precision. 

From the 'complex' approach, the width of homoclinic zone is approximated 

by 

,y+{fi)-iy-{fi) « Z,{^^) (83) 
where Z^{^^) = 4 ^o)IK ^ 

where z/g and z/4 are chosen near i/o. Since C_i(/i, 1^0) = 0{e^'^'^^) with (|8ip or 
([82]) we have 



3.9 'Real' versus 'Complex' 

The real approach provides a good estimation of the width of the homoclinic 
zone. More precisely, formula (|69p gives an estimation of the width with a 
relative error of the same order, see (j70p . However, this approach requires 
to compute the splitting determinant with the same relative error. This task 
becomes more and more delicate as the main parameter approaches 0. The 
complex approach requires less precision for the computation of the splitting 
determinant (and therefore can be computed much faster) as d is chosen larger. 
However, the estimation of the width is obtained with less precision. 

In the case of the Bogdanov map, we use similar notations: b is the slave 
parameter and a is the main parameter. The first harmonic computed with (I47|) 
is denoted by R_i(a, b), where b is the analogue of D in the case of the Quadratic 
map. Simlarly, C_i(a, 69) stands for the first harmonic computed with (j8ip or 
((82|) where 60 is the analogue of vq in the case of the Quadratic map. For 
illustration, we compute the first harmonic and the width of the homoclinic 
zone using both approaches for the Bogdanov map (7 = 8), see Figure 2. We 
easily verify that 

logio(^C_i(a,6o) -R-i(aJ)) « 2 logio(|R|-i(a J)) + logio(e'"'), 

which follows from (iTOl) and ((84| . Furthermore, we also verify that 

\C^,{a,bo)\e^-'{b',~b',) \ 

« logio(6+ -b-)+ logio |C_i(a,6o)e'"'|, 
where 63, 64 are the analogues of ^3, z/4 respectively. 



logio \Zr{a) ~ Zc(a)| 



log 



10 



26 



Example: We consider the Bogdanov map when a « 7 * 10~^. Using the real 
approach, we have logiQ(&"'" — 6~) « —1000, see Figure 2. With this approach, 
we compute Qah with a relative error of order 10"^"°", which is already a 
quite delicate task. However, from the complex approach, we can (for instance) 
choose S in such a way that e^""^ « 10^°°, see Figure 2. This way, for values of 
t e [id, ii5 + 1] , we have 

logio(0a,hoW) ~ log(Ci(a,6o)e2-*) « -300. 

Therefore, computing Ci(a,5o) with (j82p requires the computation of the split- 
ting determinant with a relative error or order 10"'^*'°. Moreover, we just need 
to find a first value of 6 = 6o such that 

logiosup||r„(t) -r,(t)|l « -300. 

teio 

However, instead of having a relative error for the width of order lO^^^^" as in 
real approach case, we obtain an estimation of the width with a relative error 
of order lO'^™. 

Now that we can compute the width of the homoclinic zone, we do so for h 
(several hundred) values of fi^^"^ and establish the set 

'H = {{^i:^ ,\og{5i)), 5i^v+{^i)-v~{^i), c< ^ii<d, i = l,...,n}. (85) 

In what follows, we describe how from the ansatz ([8]) we extract the correspond- 
ing coefficients. 

3.10 Extracting the coefficients 

Recall that the ansatz we shall consider takes the form (fTTj) where the /„'s 
satisfy ([T0| . From the set H defined by ([85|) we construct the following matrices 

^ = {^i,j)i=0,...,n-l, j=l,...,n, ^ij" = fiil^^/'^)- 

In the case of the Bogdanov map, the set of normalised data is defined in ([T5|) . 
that is the ^^^'s above are replaced by a^/^. 

Let 

a = («!, . . . , an) = A^"^ • w, 
where w = (log(5i, . . . , logSn)- Observe that 

n-l 

5Z ^iMfJ-Y^) = logf^j, Vj = 1, . . . ,ri, 

i=0 

that is the coefficients a^'s have been constructed in such a way that the map 

n 

0{"> : (0,eo)^K, (x) =^a j, (a;) (86) 

i=0 
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Figure 2: Above: (I)-Graph of log^o \e^'^^\ against the parameter a. (II)- 
Graph of log;^Q |Ci(a, 6o)| against the parameter a which essentially coincides 
with the graph of log^Q |Ri(a,6)| against a. (III)-The corresponding error i.e., 
logiQ |Ci(a, 60) — ^i{a.,b\ against a. Below: (IV)-Computation of the magni- 
tude of the homocUnic zone with the real approach (i.e., logio(-Zr(a)) against a. 
(V)-The graph of log^g \2r{a) — 2c{a)\ against a. 
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interpolates the set of data Ti. 



To illustrate our techniques, Table 1 indicates the first coefficients of the 
interpolation {h « 100) in the case of the Bogdanov map (left, 7 = 3) and in 
the case of the Quadratic map (right, 7 = —3). In the case of the Henon map, 
replacing the ansatz ([8]) by p6)) . we obtain the coefficients indicated in Table 2. 

Redoing the above interpolation for different values of 7 reveals that the first 
non linear terms in the expansion satisfies 

in the case of the Quadratic map, and 

67^ ^ 



i3i(7) = -(^yj , (88) 

in the case of the Bogdanov map. These equalities are verified with a large 
precision. More precisely, we show that ([57)1 and (|88p are verified up to the 
same number of correct digits as in (j95p when checking the extrapolation to 
zero, see section 4.3 for more details. 



4 Validation of numerical method 

To test the validity of our result, we propose three tests. To begin with, we 
test the validity of the ansatz. In what follows the experiments are presented 
in the cases of the Bogdanov map and the Henon map, but the same test can 
be applied in the case of the Quadratic map hereby confirming formula ([1] 



4.1 Extrapolability 

We claim that the ansatz ([8]) is appropriate for an asymptotic expansion of the 
width if the following criterion is satisfied. 

Assume a function G : (0,eo) ^ possesses the following asymptotics at 

00 

G{x) X ^aifi{x) 

i=0 

where {fi{x), i £ N} is the asymptotic sequence defined in (fTU)) . Define 

3fc+3 

G{3'=+3}(x)= ^ a,/,(x). (89) 

i=0 

We have 

\G{x) - G{3^-+3>(x)| = x2'=+3(^a3fc+4 +ei(a:)), (90) 
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where ei(x) ~ 0{x). From ([50]) we get 

log|G(a;)-G{3''^+^>(x)| = log |a3fc+4| + (2fc + 3) logx + log(^l + e(x) 

= log|a3fe+4| + (2fc + 3)logj; + £2(a;) (91) 

where £2(2;) = Oda;]). This impHes that the quantity log|G(a;) - G^^''+^^x)\ is 
approximatively Hnear in log a;. This must be satified for values of x outside the 
data set used for interpolation. 

Now we apply this criterion to the Bogdanov family. Recall that a is the 
slave parameter and b is the main parameter. Take an interval [c', d'] where 
c < c' < d' < d and consider the interpolation of the set Ti for values of a in 
[c', d']. In other words wc consider the set 

n' = {{a^^^,logS{aj)en\c' <a<d'} 

that consists of 3fc + 4 different values and construct the corresponding set of 
coefhcients {ai}i=o,...,3/c+3 as described in section 3. We plot the set 

£e,d = {(log(a),log|G{^''^+^>(VH) - ib+{a) - r(a))|), c < a < d} (92) 

in Figure 3: n = 140, fc = 36, c = 3.5 * 10"^^ d = 9.4 * 10"^ xhe bold line 
shows the interval [c',o?']. From and ([50)1 we must get 

log\G^^''+^^Va) ~ (6+ (a) - b" {a))\ w 37 log a + C. = 741og\/^+ C, 

where C is a constant. In Figure 3, the set looks like a straight line with a 
slope ~ 75, which indicates that the ansatz ([8]) satisfies the above criterion. 

In the Henon case, we interpolate the set of data with the polynomial 
ansatz ([9]) and the normalised width takes the form 



^^(^l^^xyi,(i-a)V4 

K{l-a,0) ^ ' 



We test the polynomial expansion the same way we test the Dulac expansion 
for the Bogdanov. More precisely, writing 

k~l 00 

G{^-i}(a;) ^^A,x\ G{x) -^A.x' 

we have 

log|G(x)-G{^-i>(:r)| =log|4| + fclogx + 0(x) (93) 
and replacing a; by (a — 1)^/^ in leads to 

log |G((5- 1)1/4) _GW((a- 1)1/4)1 ^ iog|4| + ^iog(5-l)) 

+ o((a- 1)1/4). 
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-100 



Figure 3: Plot of the set Cc.d (l92)) for the Bogdanov map in the case 7 = 3, 
c = 3.5 * 10-^ d = 9.4 * 10-3, n = 140, k = 36. 



is plotted (with k = 60, c = 1.69 * 10"^", d = 1.125 * 10"'^) in Figure 4 and 
mimics a straight line of slope « 15 = 60/4, meaning that the polynomial ansatz 
satisfies the above criteria. 

The second experiment consists of checking the stability of our interpolation 
when changing (randomly) the data TL. 

4.2 Checking numerical stability 

In this section, our interest is with the precision of our data for the normalised 
width of the homoclinic zone that is required in order to produce reliable results 
for the coefficients. The result of our test is presented in the case of the Bogdanov 
map, i.e., we test the asymptotic expansion (|15|) . In order to simulate round-off 
errors, we modify the data in the iV-th digit by adding a random perturbation of 
order 10"^ to every value of the normalised width and recompute the coefficients 
of the asymptotic expansion using the procedure described in section 3. We 
repeat the experiment for several values of N . Figure 5 concerns the coefficients 
All in pS]) : for each value of TV, we recompute the corresponding coefficient 
(denoted by ^{f^'^) after adding a random perturbation of order 10"^. Figure 



The set 




(94) 



31 



-260 



-2 8 



-300 



-320 



-20 



-19 



-18 



-17 



Figure 4: Plot of the set Cc,d O for the Henon map, c = 1.69 * 10"^", 
d = 1.125 * 10"^ fc = 60, n = 140. 

5 clearly indicates, that the precision of the computation decrease linearly with 
respect to N and if the number of correct digits in the data is less than N = 
170, the corresponding coefficient cannot be computed correctly. However, if 
N = 200 the corresponding coefficient is computed with 30 correct digits. 

4.3 Extrapolation to zero 

As announced in Section 2, for each family we are able to define the splitting 
constant associated with the 'unperturbed map', see [20| for more details. In 
what follows, our discussion concerns the Bogdanov family. The splitting con- 
stant is denoted by 0(7). Using formula we have 



Since we can independently compute the invariant ©(7) with a very high pre- 
cision, we can easily check the validity of our computation for the first term of 
the asymptotic expansion. The following table indicates, for different value of 
7 the values of 0(7) (left) computed with 20 correct digits. For each value of 
7, we observe that (j95l) holds and we indicate the relative error represented by 
— logio l(®(7) ~ exp(^o(7)))/0(7)| in the right column. 



exp(Ao(7)) 



0(7). 



(95) 
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-20; 
-40- 



160 17-0 ■■■.-...,.18 190 200 210 



1^11 -Ai 



Figure 5: The relative error logj^Q | — — — | plotted against N for the 

Bogdanov map "with 7 = 3. 



7 


0(7) 


-logio 1(0(7) -exp(Ao(7)))/0(7)l 


-2 


0.28524883190581352 


65.23 





2.47442559355325105* 10** 


90.01 


3 


4.05522622851113044* 10^^ 


62.04 


6 


2.70980378082897208* 10*' 


60.03 


7 


3.09943158275750458* 10^* 


59.6 


9 


5.18377311752952789* 10^^ 


55,6 



Table 3: The value of 0(7) for different values of 7. We clearly observe that 
the splitting constant coincides "with the first term in (jlSp 
up to the first 50 digits at least. 
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